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Abstract: 

We provide a method for the all order computation of small x contributions at the 
leading logarithmic level to cross-sections which are differential in rapidity. The 
method is based on a generalization to rapidity distributions of the high energy (or 
kx) factorization theorem hitherto proven for inclusive cross-sections. We apply the 
method to Higgs production in gluon-gluon fusion, both with finite top mass and in 
the infinite mass limit: in both cases, we determine all-order resummed expressions, 
as well as explicit expressions for the leading small x terms up to NNLO. We use 
our result to construct an explicit approximate analytic expression of the finite-mass 
NLO rapidity distribution and an estimate of finite-mass corrections at NNLO. 
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1 Introduction 



Techniques for the computation of small x logarithmically enhanced contributions 
to hard QCD cross-sections to all orders have been available for some time, and 
have been successfully applied to a variety of hard processes: heavy quark photo- 
and lepto-production pQ, Deep-Inelastic Scattering [2], and more recently hadropro- 
duction processes, including heavy quarks J3J. Higgs (both in the infinite top mass 
limit [I] and for finite m top Ej)> Drell-Yan [7] and prompt-photon [8]. In each 
of these cases, the leading contribution to the hard cross-section in the x — > limit 
has been determined to all orders in the strong coupling a s (Q 2 ), where Q 2 is the 
hard scale of the process (typically related to the invariant mass of the final state 
for hadroproduction processes), and x is the dimensionless ratio of this scale to the 
available (partonic) center-of-mass energy x = Q 2 /s. 

The interest in these results is twofold. First, they provide information on un- 
known higher-order contributions (or cross-checks on known ones): for example, the 
NNLO leading small x contribution to Higgs production in the finite m top case [5] 
was used in Ref. [9] to construct an approximate expression for the yet unknown 
full NNLO term by combining it with an expansion in l/m top [El HI]- Second, they 
can be used to construct all order resummed results, which may be necessary in 
the kinematic region where a s In - ~ 1 so that the conventional perturbative expan- 
sion becomes unstable. For many processes mentioned above, such as Deep- Inelastic 
Scattering or Drell-Yan, the Born cross-section is a quark-induced 0(oP s ) electroweak 
process (parton model), so that the leading-log resummation of hard cross-sections 
is only consistent (i.e. resummed results are only factorization-scheme invariant) if 
combined with a next-to-leading log resummation of evolution equations. Such a 
consistent resummation is possible because not only leading [12] but also next-to- 
leading [13] In i contributions to evolution equations are known since some time, 
as well as methods to use resummed evolution equations in a way compatible with 
physical constraints [HJ [T5] , and to combine solutions to resummed evolution equa- 
tions with resummed hard partonic cross-sections [HJ ELT] . 

The usefulness of these results for phenomenology has been hampered so far by 
the fact that they, and the computational techniques on which they are based, are 
only available for fully inclusive cross-sections. This is likely to become a very serious 
limitation at the LHC where, for instance, the LHCb experiment will measure Drell- 
Yan rapidity distributions at very small invariant masses (i.e. at very small x) but 
in a limited rapidity range [IB]. It is the purpose of this paper to overcome this 
limitation by extending the resummation formalism to rapidity distributions. 

Available calculations at small x are based on a suitable formulation of QCD 
factorization which is consistent with all-order leading log x resummation, called 
high energy or factorization [TJ [T9]. Generalization of this factorization to the 
less inclusive case is a prerequisite for resummation, and it will be accomplished 
here. We will proceed as follows. First, in Sect. [2] we will rederive the standard high 
energy factorization for inclusive cross-sections of Refs. [U [19]. This derivation is 



3 



completely equivalent to that of Refs. [TJ [19] , but it differs from it because the ladder 
expansion which underlies the factorization is treated in a way which is "dual" to 
it in the sense of Ref. [20J: small x logs are factored as part of standard DGLAP 
splitting functions, rather than by solving a BFKL [12] equation for the gluon Green 
function. 

The consequence of this is that the same kinematic as in the proof of collinear 
factorization (21] can be used. Within this kinematics, the dependence on transverse 
and longitudinal momentum components are kept separate from each other, and this 
renders the generalization to rapidity distributions possible. This result is derived 
in Sect. EJ where a simple formula for the computation of leading log x contributions 
to rapidity distributions to all orders is arrived at, in terms of the Fourier-Mellin 
transform of the partonic rapidity distribution computed at the lowest nontrivial 
order with incoming off-shell gluons. 

In the second part of the paper we apply this formalism to Higgs production in 
gluon-gluon fusion. First, in Sect. H]we summarize the small x behavior of rapidity 
distributions in this case, both in the limit m top — > oo and for finite top mass, and 
specifically we obtain explicit results for LLx terms up to NNLO, already known 
at NLO in the m top — > oo limit, but otherwise determined here for the first time 
in closed form. In Sect. [5] we then apply to this case the resummation formula of 
Sect. [3J by deriving a resummed expression and explicit coefficients up to NNLO. 
We check that the results up to 0(aj) of the previous section are reproduced (which 
turns out to happen in a rather nontrivial way), we construct an explicit analytic 
matched O(a^) expression in the finite m top case, and compare it to known numerical 
results [26, 27J. We also give an estimate of finite mass corrections at NNLO. 
Technical results on MS factorization and subtraction of collinear singularities both 
at the inclusive and differential level are collected in two appendices. 

2 High energy factorization of inclusive cross - 
sections 

We will rederive high energy factorization [TJ [19] in the leading logarithmic approx- 
imation (LLx) for inclusive quantities. After some introductory comments, we will 
first consider the known case of n-gluon emission at the double log (i.e. LLx-LLQ 2 ) 
level, and finally the nontrivial general LLx case. We then recall how the factor- 
ized result can be exploited for resummation, and finally we compare directly our 
approach to the original one of Refs. [JJ, [19]. The key idea of our argument is to 
exploit the factorization of the partonic cross-section in terms of a hard part and a 
ladder part, but then compute the ladder part using standard collinear factorization 
and DGLAP evolution, rather than BFKL evolution, exploiting the fact that due to 
perturbative duality (2TJJ, EH] they are in fact equivalent at the leading twist level. 
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2.1 Factor ized cross-section 

In order to simplify our derivation, we consider at first a photo- or lepto-production 
process 

V{n)+g(p)^S + X. (1) 

characterized by a hard scale Q 2 . In Eq. (CQ) g(p) denotes a gluon of momentum 
p, V(n) denotes an electroweak gauge boson of momentum n, not necessarily on- 
shell: e.g. for neutral-current Deep-Inelastic Scattering V = 7* or V = Z*, and 
S the desired final state (a quark for DIS, a dilepton for Drell-Yan and so forth). 
Generalization to hadroproduction or quark-initiated processes is straightforward 
and will be discussed In Sect. 12.61 

Our starting point is the same of Refs. [U [19], namely the observation [12] EH] 
that in an axial gauge high energy enhanced terms only come from cut diagrams 
which are at least two-gluon reducible in the t-channel (see Fig. [1]). The generic 
dimensionless partonic cross-section a can then be split in a (process-dependent) 
two-gluon irreducible part ("hard" part), and a generally two-gluon reducible part 
("ladder" part): 

f Q 2 

a= —H flu (n,p Lj pjr J jj, Rj ij, F ,a s )-L flu (p L ,p,fi Rj ij, Fj a s )[dp L }, (2) 

where the hard part includes the phase space for the observed final state S and 
the momentum conservation delta function and we have also made explicit the flux 
factor l/(2s). Here pj stands for a set of momenta parametrizing the final state of 
the hard part, which we will henceforth denote as J 7 , see Fig. |2j The measure [dpi] 
is the pz,-loop integration measure for the cut diagram Fig. [1] and will be specified 
in the following (see Eq. ([7]) below). 

Diagrams like Fig. [T] contain in general both ultraviolet and infrared divergences, 
which must be regulated. This leads to the dependence of the partonic cross-section 
on the renormalization jj,R and factorization \ip scales. The renormalization scale 
dependence can be reabsorbed in the running coupling a s (nn). Since running cou- 
pling effects are NLLx, while we are working at the LLx level, in the following we 
will omit all the dependence and consider only the dependence on fi F , which 
from now on we will call just \i. 

In order to simplify our derivation, we will assume further that the hard part 
is not just 2-gluon irreducible (2GI), but in fact 2-particle irreducible (2PI). This 
implies that we can neglect the /i^ = /z dependence in H^ u . The generalization to 
the case in which the hard part is 2GI but two-particle reducible has been worked 
out in [2] ; the extension of our formalism to it would be straightforward but we will 
not discuss it further. The functions H^ u and L^ v thus depend on the following 
variables: 

H^ = H^(n,p L ,p^,a s ) 

L» v = I/ v {p L ,p,ii,ot s ) (3) 
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Figure 1: Left: decomposition of the cut partonic cross-section in terms of two- 
gluon irreducible hard part and a reducible ladder part. Right: generalized ladder 
expansion of the ladder part. 

Since everything is on-shell and we are working in an axial gauge, we can decom- 
pose both the hard part and the ladder part in terms of conserved Lorentz structures 
time dimensionless scalar functions: 

V PI J \n-pL Q 2 J 



+Pl ( 4 _ JUL) (4 - JUL) Hll (J?l, -A ^ Q< , (4) 

PI n-pL/\P L n-p L J \n ■ p L Q 2 

PI V Pl J \P-Pl -Pi 

\Pl P-PlJ \Pl P-PlJ \P-Pl -PI J 

where Qjr stands for a set of (typically angular, see e.g. [7]) dimensionless variables 
which characterize the final state J 7 . Note that we have explicitly extracted the 
propagator factor \ jp\ from the scalar functions 

In the high energy limit Eq. ($Z§ simplifies somewhat. To see this, we work 
in the center-of-mass frame of the colliding partons and introduce the Sudakov 
parametrization 

p L = zp-k krf n = ( \f^z, j=J^- -;-k T ^) (6) 



s{l-z) VV 2 s/2s{l-z) 
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where k = (0, k x , k y , 0) is a purely transverse four-vector with k 2 = —k\ < and as 
usual s = 2p ■ n, and in the last step we have written the corresponding light-cone 
components. With this parametrization it is possible to show that the integration 
measure [dpi] for the cut diagram Fig. CD is (see e.g. [21]) 



[dpi 



dz 



2(1 -z) 



T- 



(7) 



Since the z integration goes from i to 1, we can easily identify the leading small x 
region with the z < 1 region. As for k\, it is typically bounded by the hard scale 
Q 2 . Since in the high energy regime Q 2 <C s, we have k\j s <C 1. Summarizing, at 
small x we are in the regime 



z < 1; 



<i; 



(8) 



subleading terms in z and k\js upon integration lead to power-suppressed 0(x) 
terms. 

This kinematics leads to a simplification in the arguments of the scalar functions 
H\\ i , L\\ i . We have 



Q 2 



-vl 



n-p L Q 2 



-Pl 
Wl 



Hu 



fi rp 



zs" Q 2 (l -z) 
Ts^ 2 ^ 



1 /i 2 



2' k 



a. 



l + 0(z)) 
l + 0{z)) 



^a s )(l + 0(z)). 

fvrp 



(9) 



With this in mind, and using the fact pQ that the scalar functions H± and Hu have 
the same small x behavior, and similarly Lj_ and Lu, we can write the differential 
dimensionless cross-section as 



da 



2s 

Q 2 



Q 2 o 

iljr, a, 



2s^" 



O, 

2S' Q 2 



2 - 



a, (1 + 0(2)) 



2ttLi 



l + 0(z)) 

i_ 

fvrp 



dz 



d 2 k 



T 



dz dk\ d6 
z k\ 2n 



(10) 



In Ref. P the ladder part L^pljP, fi, « s ) is computed at the LLrc level in terms of a 
gluon Green function, which in turns sums leading logs of x by iterating a BFKL [12] 
kernel. Here we observe that we may equivalently express L^p^p, £t, a s ) in terms 
of the generalized ladder expansion of Ref. [21] , and we will thus express it in terms 
of standard collinear anomalous dimensions. First, however, we discuss the hard 
part. 
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Figure 2: Graphic representation of the hard part in Eq. ([2]). Note that the hard 
part contains the momentum conservation delta function as well as the S and X 
phase space integration, but it does not contain phase space integration for gluons 
emitted along the ladder. 



2.2 Hard part 

Let us now concentrate on the process-dependent hard part. We introduce a hard 
coefficient function C defined as 

r ( Q 2 Q 2 \ r de Q 2 



(W(f H\\(^X^oc 8 )n + (11) 



o 2tt 2sz \zs' k\ 



where we have defined the projector 



= -^r (12) 



which, up to O(z) terms, selects the longitudinal part H\\ of the full H^ u . In Eq. (fTTj) 
the explicit dependence of C on Qjr is understood. 

Equation (TTTj) has a simple interpretation: C is the cross-section for the partonic 
process V(n) + g*{q) —> J 7 for an off-shell incoming gluon with momentum 

q = zp + k; q 2 = —k%. (13) 

In this interpretation, V can be thought of as the sum over the polarizations of the 
off-shell gluon. Note that 



i-2-k jn i / 



q^n u + q u n^\ _ d 



q ■ n J 2 



= IT, (14) 



i.e. the azimuthal average of performs the average over the polarizations of an 
on-shell gluon with momentum q. Equations (fTT| [T4l) then imply 



lim (c (—, %-,at s ) ) = (Jon-sheii (V(n), g{zp) F) , (15) 
k T~>° \ \zs k^ / 
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i.e. in the limit kx/Q 2 — > the azimuthal average of the hard coefficient function 
C reduces to the on-shell partonic (dimensionless) cross-section. Note that in high 
energy kinematics Eq. (J8]), unlike in standard collinear kinematics, we cannot assume 
k 2 -. <C Q 2 and thus the full dependence on must be retained in the coefficient 
function C Eq. (flTl) . However, the assumption that k\ <C s is sufficient [1] for the 
factorization of the cross-section Eq. fl2]) to hold. 

For the subsequent discussion, it is important to observe that in the limit x — > 
the hard part, being 2GI, is finite, i.e. it does not contain logarithmically enhanced 
contributions. In particular, unless the vertex of the interaction for the process 
V(n) + g*(q) ->■ T is pointlike, [TJ 



This is for instance the case for Deep-Inelastic Scattering and for heavy quark, Drell- 
Yan or Higgs production with finite m top - For Higgs production in the m top — > oo 
limit the effective gluon-gluon-Higgs interaction is pointlike, and in such a case [4J 



with Co some constant. As it is well known, and as we will see explicitly in Sect. H] 
the behavior Eq. (|T7|) leads to double logarithmic small x singularities. 

2.3 Ladder part: generalities and double log approximation 

We now turn to the computation of the ladder part. Unlike the hard part, which 
is 2PI, the ladder part may contain collinear singularities. In general, the collinear 
singularities of the partonic cross-section a, Eq. (j2J), originate from those of the 
ladder part, and from the collinear region of the integration of the momentum pi 
which connects the hard and ladder part. Collinear singularities can be subtracted 
from the partonic cross-section and factored in the parton distributions, to which 
the hard cross-section is connected by the lower line with momentum p, using the 
iterative procedure of Ref . [21] , based on a generalized ladder expansion. As already 
mentioned, after regularization and factorization of the collinear singularities the 
ladder part is a function of -J , where \x is the (factorization) scale introduced in the 
course of regularization. 

In the generalized ladder expansion, L\\ is written as a sum of terms each of which, 

after regularization and factorization, leads to a contribution proportional to ln n ^ 
to the cross-section, with n ranging from one to infinity. The n-th order contribution 
comes from a ladder diagram obtained iterating n times a kernel K{pi,pi-±, a s ) 
with % — 1,2, ... ,n. To simplify our notation, from now on we will make the 
identifications p n = pi, po = p, see Fig. [TJ Subsequent iterations of the kernel are 
connected by a pair of t-channel lines with momenta p n (see Fig. [T|) and ordered 




(16) 




(17) 
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transverse momenta k% 1 -C k\ 2 >C • • • >C k% n 



k\. Because we are interested in 



the high energy limit, reducible lines are all gluon lines: quark contributions only 
appear at NLLx and subsequent orders While we refer to Ref. [2T] for a detailed 
discussion, the way the factorization of collinear singularities is performed in the 
generalized ladder expansion is summarized in Appendix El 

We will now show how to compute the ladder part using the generalized ladder 
expansion in the double logarithmic approximation, i.e. in the very simple case in 
which the result is determined simultaneously at the LLx and LLQ 2 , and then turn 
in Sect. 12.41 to the more general LLx case which we are interested in. At LLx-LLQ 2 , 
the kernel K(pi,pi-i, //, a s ) must be linear in a s because the highest power of ln n 
in the cross-section must coincide with the order in a,. We define thus 



K l 



E 

}\trp 



,a< 



EE K 



LLQ 2 / I* 



}\rrp 



with K l linear in a s . 

Consider first the case in which the kernel K l is only inserted once. The corre- 
sponding (integrated) contribution to the cross-section Eq. (fit)]) is given by 



s ' Q 2 



o. 



1 dz 

z 



Q I -* ^ Q 

krp \ zs kn 



Q 2 Q 2 



T 



(19) 



where the bar over a indicates that collinear singularities have still to be subtracted 
from the cross-section, so this expression only makes sense at a regularized level. It 
is convenient to introduce the dimensionless variables 



X EE 



Q 2 



k'jn 

Q 2 ' 



(20) 



in terms of which Eq. ffl~9l) becomes 



Q 2 



a | x, t^,« ; 



L z J £ 

Note that in terms of these variables 



27TK 1 



Q 2 i 



a, 



(21) 



H 2 /i 2 



Q 2 £, krp 



2 • 



(22) 



i.e. the combination Q 2 £,/[i 2 is Q 2 independent, and a In™ % term shows up as 
In" S£ 

The z convolution in Eq. (|21| is turned into an ordinary product by Mellin 
transformation: 



1 rc+ico ,nr 



f(N) ee / dxx^fix); f(x) 

'0 



2m 



(23) 
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(note that by slight abuse of notation we denote with the same letter both the 
function and its transform). We get 



a 1 



Q 2 i 



(24) 



Since in the limit £ — > (i.e. k\ — > ) the coefficient function C does not vanish 
(see Eq. ffl5l) ). Eq. (|24p is collinear divergent. We use dimensional regularization 
and let d — 4 — 2e. It can be shown [2] that in d dimensions V^ u is still given by 
Eq. (fT2~j) . hence Eq. (Til)]) remains true also in d = 4 — 2e dimensions. 
The regularized expression is explicitly given b)Q: 



a 1 



d£ 



C(N,£,a s ;e) 



N \W 



a 



s ,N,e 



(25) 



where the dependence of K 1 on N is 0(e). In Eq. ( 125]) the limit £ — >■ oo is never 
reached, because k\ is bounded by the kinematics (typically k\ < Q 2 ). Note that 
we have used the assumption that the hard part is 2PI, which ensures that C is free 
of collinear singularities [211 [30]. The generalization to processes for which H is 2GI 
but not 2PI (such as Deep-Inelastic Scattering and Drell-Yan) requires also treating 
the collinear singularities in H [21 [7]. Note that the dependence on /i 2 enters only 
through the combination (fi 2 y<y s and it is thus fixed by dimensional analysis. 
Using the expansion 



i=0 



ln^ 



-e 



2! 



in Eq. (J25J) we get 



'g 2 



" — °"on-shell x 

e 



-K 1 (a. 
N K s 



+ finite, 



(26) 



(27) 



where a on _ s h e ii is the on-shell cross-section for the process V(n) + g(zp) —> J 7 , see 
Eq. f|T5l) . Equation fT2~T|) shows that 2nK 1 (a s )/N is the residue of the 0(a s ) collinear 
pole, i.e. the leading-order anomalous dimension. Because all the computation 
has been performed at the LLx level, only the leading small x contribution to the 
anomalous dimension is significant, so 



2tt 
~N 



K l (a s ) = a s 7o(iV); a s 



a„N r 



71 



7o (JV) = l/N. 

Having determined K 1 , we can substitute it in Eq. ( 125|) : 



a 



N, — , a s ; e ) = / (iV, £, a s ; ej a, 1 



Q 2 







lo(N), 



(28) 



(29) 



1 Henceforth, to simplify notation we omit terms like (47r) e /r(l — e), which in the MS scheme 
are subtracted anyway. Full details on the calculation are given in Appendix El 
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where we used the fact that the LO gluon-gluon anomalous dimension is e inde- 
pendent. Equation (129]) is the expression for the regularized (collinear singular) 
cross-section V(n) + g(p) — > J 7 + g in the small x limit. 

We are now ready to compute the ladder expansion of L at LLx, LLQ 2 . The 
regularized contribution to the cross-section when the kernel is iterated n times is 



a' 



X 



N. 



a. 



-,ai s ;e 

,2 



O.s 



Q 2 



lo(N) 



Q 2 

Sn-1 



7o (iV) 



X ... X 



Sn 
6 ( 

a. 



l+e 



C(iV,£ n ,a s ;e) x 



Q 2 



7o(iV) 



tl+e • 
SI 



(30) 



Had we used a collinear cutoff \ip instead of dimensional regularization, the or- 
dered integrations in Eq. (|30|) would give rise to the expected \ia n Q 2 /fj, 2 F term. In 



dimensional regularization, Eq. (1301) contains collinear e poles. Their MS factoriza- 
tion using the generalized ladder expansion of Ref. |21J is performed by requiring 
Eq. f )30|) to be finite after each £i integration. This is done subtracting iteratively 
in such a way that Eq. f l30|) is finite also at the integrand level; details are given in 
Appendix [A] This then leaves a single n-th order e pole in the cross-section, which 
can be subtracted using standard MS. 

The result after iterative subtraction of the first n 



1 singularities is 



a"[N,E as]e 



2 \ e 



X 



£1+6 



C(iV,£ n ,a 8 ;e) 



(n - 1)! e™" 1 



« s7 o(iV) 1 - 



Q 2 ^ 



n-l 



(31) 



which contains a single l/e n pole, as it ought to in the MS scheme. 

The full LLx-LLQ 2 result is found by adding up all contributions Eq. ( 13TT) . It 
turns out that (as well known from Ref. pQ for the general LLx case) the sum is 
actually finite, i.e. the exponentiation of the leading collinear singularities to all 
orders leads effectively to a cross-section evaluated with off-shell gluons, which is 
thus finite. Indeed, we get 



a 



£ 

n=l 



Q 2 



2 \ e 



lo(N) 



d£ n(AT . « S 70(V)I 

— C*(iV,£,a s ;e)e 



(32) 



which is finite in the e — )■ limit. Indeed 



lim — — 

e-s-0 £!+ e 



cxp 



a s 7o(A0- 



1 



Q 2 £ 



cxp 



Q2 

a s 7o(iV) In — 



(33) 



and since a s jo(N) > also the last collinear divergence in Eq. (|32|) is regulated. 

We can then safely go back to four dimensions. For completeness we reintroduce 
the argument of the running coupling a s , although running coupling is a NLLx 
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effect. We obtain 

a (n, ^,a s (fi 2 )) = a s (^ 2 ) lo (N) x 

x / d^C^^-'e^^^^CiN^^s^ 2 )), (34) 
Jo 

or, using renormalization group invariance to set /i 2 = Q 2 , thereby absorbing the 
explicit n 2 /Q 2 dependence into the running coupling a s (Q 2 ), 

POO 

a (N, a s (Q 2 )) = a s (Q 2 ) lo (N) / ^(^boW-i^ ^ ^ a-(g2) j ? (35) 

Jo 

where we have made explicit the fact that the partonic cross-section no longer de- 
pends on n 2 IQ 2 . At each fixed order in a s Eq. ( 1351) reproduces the small N limit of 
the coefficient function, as can be checked by integrating by parts and expanding in 
powers of a s . 

We see that there are two sources for the iV-dependence of the cross-section a: 
the a s 7o(iV) terms and the explicit iV-dependence of C . Let us first assume that 
C(N, £, a s ) is regular in N = (i.e. C(x, £, a s ) — > for x — > 0, see the discussion at 
the end of Sec. 

C(N, £, a s ) = C(0, e, a s ) + NC'(0, £, a.) + ... (36) 
In this case Eq. ( 1351) can be written as (see Eq. fT25]) ) 

a(N,a s ) = ^JdZ&- 1 C(0,Z,a t )(l + O(N)) 

oo _ ■ 

= «^c^- LLQ2 (|) (1 + 0(N)), (37) 

i=0 

where the power k is pro cess- dependent and is encoded in C (e.g. k = 1 for 
Drell-Yan and DIS, k = 2 for Higgs). The expansion in a s /N leads to the de- 
sired a s In - expansion in x space. From Eq. (!37|) it is then clear that only the 
iV-independent term C(0, £,a s ) is relevant for a LLrc calculation. Things are differ- 
ent when C(x, ^, a s ) — >■ C 7^ for x — > 0. This is for instance the case for the Higgs 
coefficient function in the infinite m top approximation jl]. In this case C(N,^,a s ) 
is not regular in iV = 0, and the full iV-dependence must be kept. The singular 
behavior of the coefficient function in iV = leads to double log behavior in x space, 
as we shall see explicitly in Sec. El 

Eq. (135 p reproduces the LLx result of Ref. [I] in the LLQ 2 case: the LLa; coeffi- 
cient function is found evaluating the partonic cross-section with incoming off-shell 
gluon, performing a Mellin transform with respect to the gluon virtuality, and then 
letting M = 7, where M is the Mellin variable which is conjugate to the gluon virtu- 
ality, and 7 is the relevant anomalous dimension, which in the LLQ 2 case coincides 
with the small x limit of the leading order anomalous dimension 70 . 
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2.4 Ladder part: LLx case 



The argument presented in Sect. 12.41 actually also works in the general LLx case, 
at the price of some technical complication. In this case the kernel K contains 
contributions to all orders in the strong coupling a s . The way the ladder expansion 
is constructed when the kernel K is defined at a generic order in a s was explained 
in Ref. [21], and it is that which is used to compute anomalous dimensions at any 
perturbative order. In a nutshell, the ladder diagram of Fig. [1] is still constructed 
with strongly ordered transverse momenta in such a way that all large logs of Q 2 
arise due to integration over the momenta pi. Unordered contributions, as well as 
finite parts due to iterations of lower-order kernel after subtraction of collinear poles, 
are included as higher order contributions to the kernel K, which is now no longer 
2PI, even though it is still free of collinear singularities. The total cross-section is 
still the sum of contributions with n iterations of the kernel, with n ranging from 
one to infinity. 

With a generic kernel, the steps leading to Eq. (1251) remain unchanged, and lead 

to 



a 



o 



Q 2 

%C(N,£,a a ;e)K(N, 1 '' 



Q 2 i 



a.,; e 



(3f 



Because K is free of collinear singularities, the analogue of Eq. (|27|) also still holds: 

,2 x ! 



-Oon-sheii-fr' (N,l,a s )+ finite, 



(39) 



which again implies that K is just the standard anomalous dimension: 



K N 



2 \ e 



«,;e =1\N. 



2 \ e 



Q 2 i 



, OLe . 6 



(40) 



At the LLx level at which we are working 7 is the (rf-dimensional) LLx anomalous 
dimension. 

However, unlike in the LLQ 2 case, the anomalous dimension 7 now explicitly 



depends on e; this affects the way the MS subtraction is performed. Indeed Eq. ( 130 
now becomes 



x 



7 (N, 



Q 2 
/i 2 

Q 2 £n-1 



7 In, 



tl+e 
Sn-1 



2 \ e 



Q 2 in 



X ... X 



, a.,; e 



7 N, 



tl+e 
Sn 

2 \ e 



C(N,£ n ,a s ;e) x 



Q 2 i 



l+e> 



(41) 
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which, after iterative subtraction of the first n — 1 singularities becomes 



Q 2 



2 \ le ~ 



(n-l)\e n - 1 
where we have introduced the expansion 



r/j(N, a s ;e) 
Ji(N, a s ;0) 



n-l 



(42) 



T In, 



2 \ e 



Q 2 e 



,a s ;e 



Qf.; e 



i=0 



2 \ le 



Q 2 H 



(43) 



If 7 is given at the fc-th-perturbative order the sum goes from to k. 

Again, the total cross-section is rendered finite by exponentiation of the collinear 
poles: 



n=l 



x exp 



2 \ 6 



Q 2 



ls-^ ^jN, q a ;Q) 
e 4^ i 



C(N,£,a s ;e) x 



/i 2 V e ji(N, a s ;e) 



However, we must keep track of the e dependence in 7*. Expanding 

7<(JV, a s ; e) = 7<(JV, a.) + e^(N, a.) + e%(N, a s ) + ... 



(44) 



(45) 



we get 



lim — — exp 



l^ ^(N,a s ;0) 



2 \ 16 z, 



7i(iV, a s ;e) 



y(N,a s )-l 



exp 



Q 2 



TZ(N, a s ) = exp 



j(N,a s )\n ., 



Ji(N, a. 



K(N, a s 



(46) 



In comparison to the LLQ 2 result Eq. (13"21 the process independent but scheme 
dependent 71 factor has now appeared. Note that, because we are working in the 
LLx approximation, a s does not run. At NLLx and beyond the running of the 
coupling generates further contributions to 1Z, that depend on higher order terms 
in the expansion (l4"5l) in powers of e [3TJ [32]. The LLx cross-section is finally 



a(N,a s (Q 2 )) =7 



a s (Q 2 
N 



C(N,^a s (Q 2 ))K 



a s (Q 2 
N 



(47) 
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where without loss of generality we have set /i 2 = Q 2 (thus omitting the explicit 
dependence on fi 2 /Q 2 ) and we have noticed that the LLx anomalous dimension 
which enters Eq. ( )47l) only depends on a s and N through the combination As 
in the LLQ 2 case, if C (N, £, a s (Q 2 )) is regular in iV = 0, then only C (0, f , a s (Q 2 )) 
is relevant at LLx (see discussion around Eq. ()36l)). 

However, there is a further source of scheme dependence of our result. Indeed, 
we have seen that after exponentiation the last collinear singularity disappears; it 
is however present if the result is expanded out and determined order by order in 
perturbation theory. It is thus not obvious that the MS subtraction commutes with 
exponentiation, and therefore the result Eq. fj^BI) might be defined in a scheme which 
differs from MS by a further LLx scheme change, which amounts to a redefinition 
of parton distributions by a factor N{a s /N) [22l [23| [2^]. 

Hence, our result has the general form 

(48) 

with 

The factor M (and thus R) which fixes the factorization scheme has been determined 
for MS through a computation of the LLx gluon Green function [T9] . 

2.5 Cross-section resummation 

The LLx cross-section Eq. ( HHj) can be written in terms of the M-Mellin transform 
of the coefficient function C 

POO 

h(N, M, a s )=M d^-'C (N, £, a.) R m (M) , (50) 
supplemented by the condition 

M = 7(|), (51) 

namely 

a(N,a t ) = h^N, 7 (^j,a a y (52) 

Explicit resummed expressions can be obtained from Eq. ( 152]) by substituting in it 
the well-known [121 12H] expression for the LLx anomalous dimension 

-> GO = I + 2C(3 > (I) 4 + 2 « 5 > (I) 6 + 12C(3 > 2 (I) 7 + - (53 > 

and for R the appropriate to the factorization scheme. 
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The result (I50H51P coincides with that of Ref . [2] , where it was derived by solving 
the small x evolution equation for the gluon Green function. Indeed, our derivation 
parallels that of Ref. [2], the difference being that we have determined the ladder 
part using the generalized ladder expansion of Ref. [21], while in Ref. [2] it was 
computed as the gluon Green function. These coincide because of DGLAP-BFKL 
duality [20]. 

Explicitly, our Eq. f jTQj) in Ref. [2] is written in the form 



°= l d i d -j§C(-,^,a s )g(z,k* T ), (54) 



where Q is a gluon Green function, which satisfies a LLx (BFKL) evolution equation 
and C is defined as in Eq. (jTTI) . Taking an M-Mellin transform Eq. (150|) both of the 
coefficient function C and the Green function Q, and noting that in M space the 
Green function has the form 

G(N, M) = r( f° Mn [1 + 0(M - rta./N))} , (55) 
M - 'j{a s /N) 



one ends up with the expression 



c+i °° dM ( /j 2 ^ ~ M 



*(N, a.) = / — I ^ ) h(N , M)Q{N, M) ( W ) 



of the resummed cross-section, which is immediately seen to coincide with Eq. 
(with /i 2 = Q 2 ), by performing the M integral and identifying the R function with 

R(a s /N) = r( 7 (a./JV)). (57) 

The form (1541) of the resummed result has the structure of a high energy fac- 
torization theorem. It is important however to note that a on the left-hand side of 
Eq. (|54"|) is a partonic cross-section. Namely, the basic insight of Ref. p] is that, 
after factoring the hadronic cross-section into parton distributions convoluted with 
a partonic cross-section, at the LLx level the partonic cross-section can in turn be 
itself factorized into a hard part and a ladder part. The argument presented in this 
section has pushed one step further the analogy between the factorization of the 
hadronic process and the subsequent factorization of the partonic cross-section, by 
using the standard collinear factorization technique of Ref. [21], albeit at the LLx 
level, to compute the ladder part of this second factorization. 



2.6 Hadroproduction and quark-initiated processes 

We now discuss briefly how our results can be generalized to the case of hadropro- 
duction processes. This differs from the case discussed so far because there are now 
two incoming partons. Consider the process gin) + g(p) — > <S + X. We have now 
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two different transverse momenta kr, one for each leg. In principle, the interfer- 
ence between the two legs could spoil the factorized result Eq. (j2J); however, these 
contributions do not appear at LLx [lj, so the result in this case is simply found 
by applying the generalized ladder expansion to both legs. In particular, Eq. (j2J) 
generalizes to 

f Q 2 

where again we have assumed that H is 2PI, fi is the factorization scale and we have 
omitted any dependence on the renormalization scale. We have assigned bar index 
to contributions from the upper leg. In Eq. ( 1551) n L is defined through the Sudakov 
parametrization (see Eq. (jS]) for comparison) 

k 2 

n L = zn-k T p, (59) 

s(l-z) 

with k = (0, k x , k y ) and k 2 = —k\ < 0. In this parametrization we have (see Eq. ([7])) 

[dn L ] = ^—f^ = ^—fTdQ (60) 

In the high energy regime we have z<Cl, (see (j5J). Given the factorized 

form Eq. ( 155|) . we can define a hadroproduction coefficient function C: 

c(-,Z,£,a s ) = ^^-V^V-, p H^(n L ,p L ,p T ,a s ), (61) 

\ZZ / Jq Z7T Z7T ZSZZ 

where again we have omitted the Qjr dependence and 



L.2 h^b C 



Eq. (16T1) should be compared to its photo- or lepto-production counterpart Eq. (1TT1) . 
In particular, we now have (see Eq. ( 1T5|) ) 

lim (C ( — , f , f , a s ) ) _ = a on _ she ii {g(zn),g(zp) ->• J 7 ) . (63) 



Starting from the hard coefficient function Eq. ( 16T1) . we can repeat the same 
steps of the previous subsections, and arrive at the hadroproduction formula (with 

^ = Q 2 ) 



a 



x f & 7 (^) fW-Cl^.Ulfll) , (64) 



to be compared to its photo- or lepto-production counterpart Eq. (|48|) . 

Before concluding this section, we now briefly account for quark-initiated pro- 
cesses. The analysis in [12] shows that at the LLx level only gluons enter in the 
generalized ladder, with one exception: since the Pg q splitting function is singular 
at small x, the first rung of the ladder can be a quark. In other words, the initial 
quark can emit a collinear quark, thus becoming a gluon which starts radiating the 
generalized ladder Fig. [Tj, leading to a cross-section proportional to Pg q . Since at 
small x we have Pg q = ^Pg 9 , the only difference between quark- and gluon- initiated 
processes is a color-charge factor. 



3 High energy factorization of rapidity distribu- 
tions 

The formulation of high energy factorization of Sect. [2] will now be extended to 
rapidity distributions, thereby leading to a method for the determination of LLx 
contributions to all orders. The argument will closely follow that of Sect. [2J 



3.1 Kinematics and factorization 

We now concentrate from the onset on hadropro duct ion of a final state system S 
with invariant mass Q 2 : 

h 1 (P 1 ) + h 2 (P 2 )^S(p s )+X (65) 
and we study distribution in the rapidity of S 

_ 1 E s +p Sz 

where Eg, p$z are measured in the hadronic center-of-mass frame. The (dimension- 
less) hadronic rapidity distribution is written in terms of parton distributions 
and a partonic rapidity distribution -^ L : 



{x h ,Y,Q 2 J a s ) = ^ / dx x / 



dcr h 
dY 

a,b 



dX2 X 

V% h e Y J V*h e ~ Y 



dy \xxX2 Q 2 2 x 2 

where Xh = ^- and S = (Pi + P2) 2 ■ Note that we have used the well-known fact 
that in collinear factorization the hadronic and partonic center-of-mass frames are 
related by a longitudinal boost. 
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Henceforth, as in the previous section, we work in the partonic center-of-mass 
frame, where X\ = x^ — 1 (thus the partonic rapidity y is equal to the hadronic one 
Y) and the partonic rapidity distribution is given by 

do f ^ \_f |Mp / 1 Es+psA 

where the bar stands for sum (average) over final (initial) helicities and colors and 
dUf is the final particle phase space, and here and henceforth for simplicity the 
parton indices a, b have been omitted. 

Resummation at the differential level in rapidity can now be performed by re- 
peating the argument of Sect. [U but with the partonic cross-section now multiplied 
by the kinematic constraint enforced by the rapidity delta in Eq. ([6"8j) . In the fac- 
torized expression Eq. ( 158|) the dependence on the momentum p$ of the final state 
system S (see Fig. [2]) is entirely contained in the hard part. Therefore, the kinematic 
constraint only enters the definition of the coefficient function C which now becomes 

n ( - cc \ f 27rdd [** S Q 2 

x [V, u V, p H^(n L ,p L ,p T , a,)] sfy-Un Es + Ps A , ( 69 ) 

V z As — PSz J 

where z, z, £, £ were respectively defined in Eqs. (0 EH EH and we have made 
explicit the dependence on the rapidity y, which in turn is fixed in terms of the 
remaining kinematic variables by the delta in Eq. fl68|) . Note that, as in the inclusive 
case Eq. fl6T|) . in general the coefficient function C y also depends on a residual set Qjr 
of final-state variables which survive the phase space integration. This dependence 
is omitted for simplicity. 

Equations f lBTj I6T?|) show that the ladder part is unaffected by the rapidity con- 
straint. The rapidity-dependent coefficient function C y instead acquires a nontrivial 
dependence on all kinematic variables because of this constraint. However, we will 
now show that invariance under longitudinal boosts of the ladder implies that the 
dependence of C y on the longitudinal momentum fractions z, z is completely de- 
termined by its dependence on the rapidity y and the scaling variable x. But the 
argument presented in Sect. [2] shows that the effect of the ladder insertion on the 
external legs (which generates the LLx series) can be cast in terms of QCD evolu- 
tion in fcy, characterized by standard anomalous dimensions. Therefore, this implies 
that we will be able to trade the rescaling of the momenta pi and Ul which enter 
the hard part for a shift in rapidity. This will enable us to determine fully the 
LLx resummation of the rapidity distribution from the knowledge of the rapidity 
dependence of the hard part C y . 

In the remainder of this section, we will prove these results, which are essentially 
of kinematic origin, while the subtraction of collinear singularities in the ladder part 
is unchanged and can be handled as in Sect. [2J Therefore, for clarity we will simply 
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Figure 3: Single emission kinematics from both legs in terms of light-cone Sudakov 
components, see Eq. (jfJJ). 



regulate collinear singularities with a cutoff. The price to pay for this is that we 
will loose full control of the factorization scheme. The full computation in the MS 
scheme is presented in Appendix [Bj and will provides us with a determination of 
the scheme- dependent factor R. 



3.2 Ladder part and rapidity dependence 

As we have seen in Sect. |2j large logs of - come from gluon ladders attached to the 
external legs of the hard part. We now show how to compute the effect of these 
ladders on rapidity distributions. The structure of the computation is unchanged: 
the ladder is constructed by iterating n times a kernel K, which is closely related 
to the standard QCD anomalous dimension. In the case of single insertion of the 
kernel on the lower leg (say) we have 



da 1 f fi 



,2 



dy V ' Q 2 



dZ f dkrp fx krp f krp \ \ . /nn,\ 

W V I I' Q2 ' y [ X ' Z ' Q2 J ' a s) K a s) , ( 70 ) 

which only differs from Eq. (fT9|) of the inclusive case because of the replacement 
of the inclusive coefficient function C with the differential one C y , and because of 
the kinematic constraint imposed by the delta function Eq. (|68|) which relates the 
rapidity y to the other kinematic variables of the hard part. 

The kinematics of the single insertion is displayed in Fig. [31 using the notation of 
Eq. ([6]), in the LLQ 2 case in which the kernel corresponds to single-gluon emission 
as discussed in Sect. 12.31 As discussed in Sect. I2.4l the kinematics in the general LLa; 
case is the same. It is apparent that 



Vl=[\I~zA-Kt) (l + Ol^)). (71) 
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It follows that the longitudinal components of pi are related to those of the Sudakov 
base momentum p by a longitudinal boost. Therefore 

C y (x, z, kj., Q 2 , y (x, z, , a s ) = C y 0, y (x, - ^ In z, a s ^j +0 O^J , 

*(72) 

i.e., the effect of a i^-kernel insertion on the rapidity dependence of the hard part 
is a longitudinal boost, up to subleading terms in small x kinematics Eq. (JHJ). 

Eq. ([TP]) can be easily generalized to the case of n-kernel insertion, with the 
result: 

da n ( fi 2 \ f°° dk 2 1 , n _, k 2 
'x,^-,y,a s ) = -£-. — ln nl ^x 



dy V Q 2 ' ' / Jfi 2 k\ — 1)\ jj? 

— P(z n ,a s ) x ... x f —P(z 1 ,a s )x 



z 2 ...z n 



x °v (i^' qi> y ~ \ ln Z1 - - ~ \ ln * n ' as ) ' (73) 

where P is the inverse Mellin transform of the kernel K and for simplicity we have 
written y(x, k^/Q 2 ) = y. 

Equation fl73|) can be factorized by taking a Fourier-Mellin transform, defined 
(together with its inverse) as 

/*1 /*oo 

f{N,b) = / dxx"- 1 / dye»»f{x,y); 

JO J-oo 

fC+iOO TAf pOO TT 

/(*,*) = / ?le-*"Mb), (74) 

Jc-ioo 2m J-oo 2lT 

where as usual by slight abuse of notation we use the same symbol to denote the 
function and its transform. We get 



/CO 1 /* 

dy e iby / dx x N ~ l x / 
oo JO Jx 



1 — e i t ln -P(z n ,a s )x...x 

1 



z 2 ...z n 



x k 2 



z\ J ^ k A T [n - 1)! 



XC *{-^>V> V ' a ')- (75) 
Each z integral in Eq. (!75|) is a convolution, hence in Fourier-Mellin space the result 
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factorizes: 



da n 
dy 



7 [N + 



ib 



dkri 



b 2 

fj, 2 ™T 



Cy[N,-^,b, 



ln n-i _T1 
Ml 

(n- 1)!' 



(76) 

As in Sect. [2J the final result for the cross-section is obtained by adding up the 
contributions from n insertions of the kernel, with 1 < n < oo. We get 



da 



^da^ 

n=l y 



7 



dkrp(k< 



2 ^(jv+ii/a) 



-1, 



(77) 



3.3 Resummation of rapidity distributions 

Equation (177|) provides the resummation of LLx terms due to radiation from one of 
the two outer legs of the hard kernel, so a full resummation can be obtained from 
its two-leg generalization. Before discussing this, we observe that in Eq. (1771) the 
multiple poles in iV = which in the inclusive result Eq. (I46jl lead to the LLx terms 
have now become multiple poles in N = —ib/2. On top of these, the hard coefficient 
function C y may still contain poles in N = 0, but, as discussed in Sect. \2.2\ this 
only happens in the somewhat pathological case of pointlike interactions. 

In order to understand how poles in N = —ib/2 can be related to large logs of x 
we note that the value of x determines the range of the partonic rapidity: 

-In a; < y < -In-. (78) 
2 ~ y ~ 2 x K J 

It follows that as the energy increases, i.e. as x decreases, the rapidity range widens. 
This implies that high energy enhanced terms are not necessarily logarithmic at the 
level of rapidity distributions: for instance, a flat rapidity distribution leads to a 
contribution enhanced by In x in the inclusive cross-section because of this widening 
of the phase space. This can also be seen in Fourier-Mellin space, by noting that 
(see Eq. f !74p ) the Mellin-space inclusive cross-section is found letting b = in the 
Fourier-Mellin space rapidity distribution. Therefore poles at iV = in the inclusive 

result can be obtained from any term of the form (^ N "j^ such that /(0) = 0. 

We now turn to terms coming from emissions from the other leg. These are 
easily determined: the only difference in comparison to Eq. flTOj) is that the boost 
has now parameter 1/z (see Fig. [3]). The result is then identical to Eq. ( 1771) . but 
with the replacement: 

7 (ivT^) 7 (a 7 ^) ■ (79) 
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The resummation prescription is thus the following. First, the rapidity distri- 
bution for the process g*(pi) + g*{ni) — > S is determined with off-shell incoming 
gluons. The rules for computing C y are the same as presented in Sect. 12.21 for the 
inclusive case, namely, the incoming off shell momenta are pl = p + k, p\ = —k\ 
and til = n + k, n\ = —k^, and the average over the initial gluon polarizations is 
performed by the projectors V^ v Eq. f lT2|) and P^p Eq. fl62|) . The coefficient function 
C y is then defined using Eq. (1691) . 

In terms of the differential coefficient function C y the resummed rapidity distri- 
bution is given by the two leg generalization of Eq. ( 1771) : 



X 



oo 



I W~S C^r^^C, tN,*L,*L,i ia .) . (80) 



M 2 V 2 



P P 

IX rp IX rp 



This result has been obtained with a collinear cutoff /i 2 . Although the full MS 
derivation is given in Appendix |Bj here we give the final result: 

x {wtw) R ™ {wtw) r^'-'Ci (JUXM.) , (81) 

where £ = k^/Q 2 , £ = k\jQ 2 and we have chosen a s = ot s {Q 2 ) and consequently 
dropped the explicit fi 2 /Q 2 dependence. 

We can cast Eq. ( 18~TT) as an M-Mellin transform of the off-shell rapidity distri- 
bution C y (N,£,£, b, a s ) in a form similar to that Eqs. ( 1501 15T1) of the inclusive case, 
by defining a double M-Mellin transform of the differential coefficient function C y . 
Setting again /i 2 = Q 2 we get 

/»oo /*oo 

h y {N, M, M, b, a,) = MM / c^" 1 / fi^Cy^N, £, £, 6, a s )R{M)R{M), 

(82) 

with again a s = a s (Q 2 ). The resummed result is then found from Eq. (182]) through 
the identifications 



M = 7 



O.s 



AT+f 



^ = 7 ), (83) 
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namely 



% [N.i.a.m = h, (n,, , 7 . (84) 

The inclusive result (i.e. the two leg generalization of Eq. ( 1521) ) is recovered from the 
differential one Eq. (I84p by setting 6 = 0. Resummed results at the hadronic level 
are finally obtained by substituting the resummed differential coefficient function 
Eq. (I84p in the factorization formula Eq. (|67|) . 

As in the inclusive case, the factorized expression Eq. (IHTj) can be viewed as a 
factorization theorem, by rewriting it as 

However, a direct derivation of Eq. ( 18 5 p from the determination of the gluon Green 
function and iteration of the BFKL kernel from its inclusive counterpart appears 
to be less immediate than our derivation. This is because in our case the rapidity 
dependence along the ladder, which has the standard kx ordered structure, is trivial, 
while the rapidity flow along a BFKL ladder is less obvious. 



4 Higgs rapidity distribution at small x 

As a first application of the formalism developed in the previous part of the paper, 
we will consider Higgs production in gluon-gluon fusion. This process is particularly 
simple kinematically because the final state J 7 contains just one particle, the Higgs. 
In contrast, recalling from Sec. [2] that J 7 is the lowest-order final state for the given 
inclusive process with at least one incoming gluon, it is clear that for Drell-Yan [7j 
or prompt photon [8] production J 7 is a two-particle final state. 

In the Standard Model the Higgs couples to gluons via a heavy quark loop. If the 
top mass m top is much bigger than the Higgs mass one can describe the interaction 
through the effective Lagrangian 

dm = ^f^ HG^G^ (86) 

which corresponds to a pointlike interaction between the gluons and the Higgs. As it 
is well-known [33] this approximation is reliable within 10% for Higgs masses below 
1 TeV. However, at small x this approximation fails: when the partonic center- 
of-mass energy squared s is large enough (i.e. at small x) gluons can resolve the 
quark loop. Consequently, the effective theory at small x leads to spurious double 
logarithmic singularities [U H] , while the full theory has single log behavior: 

„„„ X J 5(l-x) + J2 k c k a k s \n 2k - 1 l m top -> oo 

° I S(l-x) + j: k d k a k s \n k - 1 l finite m top [ * ' } 
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(a) (b) (c) 

Figure 4: Graphs contributing to the small x Higgs rapidity distribution at NNLO. 
Here the (orange) blob stands for a HG a ^ v G a '^ v effective coupling (heavy top limit) 
or for the top quark loop (full Standard Model with finite m top ). Only relevant 
light-cone fractions (in units of a/s/2) are reported. 

The treatment of rapidity distributions for this case is particularly simple due to 
the fact that, using the small x kinematics Eq. flH]), the rapidity is entirely fixed in 
terms of the momentum fractions of incoming partons. Indeed, for single emission 
(see Fig. [3]) we have 

k 2 

y=\\^ = \\n l -^^ = \\n^ + 0( k ^)+0{z) (88) 
2 p H 2 z 2 z \ s J 

if the emission is from the upper leg and 

y = ~ In Z k% =ll nz + +0 (^\+0 (z) (89) 

1 s(l-z) 

if the emission is from the lower leg. In case of a single emission from each of the 
incoming legs (see Fig. H]-b), we have 

^=2 ln --V = 2 ln i + 2 ln r^fe' (90) 

with Q 2 = mf^ (hence x = m 2 H /s) for Higgs production. But using momentum 
conservation 

m% 1 (gi) 



zz 



in 



I + (k T + k T ) i + e + e + 2 vee cos ^ 



so £x/(zz), £x/(zz) is at most 0(1) and therefore 



y = Un Z - + 0{z,z). (92) 
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Equations (I88ti92p show that for Higgs production the cross-section at fixed y is just 
a cross-section at fixed z, z. Of course, these purely kinematic relations hold for any 
2—7-1 process. 

In this section, we compute the small x behavior of the Higgs cross-section (in the 
gg channel) up to NNLO, before turning to the all-order resummation in the next 
section. In the m top — > oo limit, we will obtain analytic closed- form results which at 
NLO may be compared to the known [25] full analytic expression. For finite m top 
we will reduce the NLO and NNLO small x coefficients to quadratures and tabulate 
them numerically as a function of r (i.e. of the Higgs mass, see Eq. ( 11 19ft ). In 
this case a direct comparison to the known NLO result [31] is not possible because 
the latter is only known in the form of the numerical integral of a fully differential 
cross-section. 



4.1 Heavy top approximation 

In the infinite m top approximation the NLO Higgs rapidity distribution is known in 
closed analytic form [25]. Here we reproduce the NLO result in the small x limit. As 
we have mentioned in Sec. [2], only ladder-type diagrams contribute to LLx. Hence 
we must compute the two diagrams Fig. [3j We start with emission from the upper 
leg. Up to power-suppressed terms in z and k^/s, the fully differential cross-section 
is 



a z 2 s 2 2tt . 
x — ^rd 



2sz 



a -d 

x 



m 



1- Z - + 



H 

x 



rn 



a s n 2t N c dz dk\ (47r) e 



7T 



(4) 



l+e 



r 1-e 



_ dz d£ (Any 



(93) 



where we have defined 

ct = G F V2a 2 j57&ir (94) 

and in the last step we have set fi 2 = Q 2 = m 2 H for simplicity. As in Sec. [2} the bar 
over a indicates that it contains an unsubtracted collinear singularities, so Eq. ( 1931 
only makes sense at a regularized level. 

Note that the hard part, contained in the first square bracket in Eq. (|93|) . reduces 
to the LO result when kj, — 0, z — 1, as it ought to according to Eq. (|15p . The 
only kx dependence of the hard part comes from the phase space delta function, i.e. 
the hard amplitude is k? independent. This makes high energy factorization almost 
trivial in this case. Indeed, to compute the rapidity distribution, we just have to 
insert in the differential cross-section Eq. (193]) the rapidity 5 function (see Eq. (|68|) ). 
which in our case is 



1 2 p~ H 



5 f y In - 

V 2 z 

2u5 (z — u) , 



(95) 
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where, following Ref. [25] @, we have defined 

u = exp(— 2y) 



(96) 



and the first equality in Eq. ( 195]) holds up to subleading terms. In terms of u the 
kinematic limits Eq. ( !7H|) on the rapidity become 



x < u < — . 

X 



(97) 



The transformation y -H- —y corresponds to u -H- 1/u symmetry. For the emission 
from the upper leg one obtains the x < u < 1 part of the rapidity distribution (as 
it can be immediately seen from the equality u = z), while emission from the lower 
leg covers the other half of the rapidity range. 

Using the rapidity S function to perform the z integration and the momentum 
5 function to perform the £ integration we immediately obtain the Higgs rapidity 
distribution clS db function of u: 



du 



(47T) 



1 da 

r~ — o"oQ^s„ 

2udy T(l 



x 



[u — X 



Me' 



(9* 



where we have now omitted the explicit dependence on a s for simplicity, as well as 
the dependence on fi 2 /Q 2 because we are working in the case /i 2 = Q 2 = m 2 H . For 
e = Eq. (I98p has an endpoint rapidity singularity in u = x. We regularize it by 
defining the plus-prescription 



du- 



/(«) 



u — x) 1+€ 



VM-/M +/M 



(u - x) 1+e 
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du- — 

(u 



X 



duf(x) 



1 



[u 



~-8(x - u)(l - x) e 



(99) 



where the distribution is in u space. Using Eq. 
traction our result reduces to 



99]) and performing the MS sub- 



leg . 



du 



(u — X) 



— S(u — x) In x 



(100) 



up to terms which vanish like powers of x. In Eq. f llOOp we have taken the limit 
e — > (after MS subtraction). Integrating Eq. ( 1 1 f) over u, the first term gives a 
vanishing contribution, while the delta function integration leads to o = <Joa s In 1/x, 
in agreement with the inclusive result [I]. 

One may be tempted to think that in the small x limit the first term on the right- 
hand side of Eq. f llOOp can be neglected in comparison to the second one. However, 
this is incorrect because the first term is phase space enhanced. To see this, we use 



2 Note that even though we follow closely the notation of Ref. [25], we depart from it in calling 
y the rapidity, whereas in Ref |25j y denotes a certain function of u and z. 
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Fourier-Mellin space. Start with the first term on the right-hand side of Eq. (11001) . 
Its Fourier transform with respect to u is 



/ ib 
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2i 1 


/ U 2 








_{u-x) + _ 





' 2m 
+ m+ 
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-X 



2 F 1 (l,l + f,2 + f,x) 



1 + 



ib 



x 2 In x 2 ln(l — x) — x 2 ip ^1 — — j . 

(101) 



Note in particular that the In 1/x term in the last line exactly cancels the Fourier 
transform of the delta S(u — x) In 1/x contribution in Eq. fllOOp . Performing the 
Mellin transform, we get 



{u - x)+_ 



1 - NH{N) + NH(N + f ) 
N(N-%) 



(iV-f) 



6\ 2 ' 



(102) 



where H(N) is the iV-th harmonic number. Turning now to the second term on the 
right-hand side of Eq. f llOOl) . its Fourier-Mellin transform is 



S(u — x) In 



x 



(iV-f) 



2 ' 



(103) 



Combining the two terms we get 



(u — x) + 



5(u — x) In x 



1 - NH(N) + NH (N + f ) 
N(N-%) 



+ 



l-NH(N)+NH(N+%) l + 0(N) 



(iV-f) 2 

(104) 



(AT- |) 2 "" N(N-f) "iV(iV-f)' 

The inverse Fourier-Mellin of the leading term in this result leads to the cross-section 



da i-ieg 
du 



[x,u) = MT 1 



iV(iV-f) 



cr a s -, 
u 



(105) 



to be compared to Eq. fllOOH . Away from the endpoint x = u, Eq. f llOOl) admits the 
small x expansion 



d<j\. 



du 



I L / x f x 

(x, u) = a a s = cr a s - 1H h - 

u — x u \ u \u 



(106) 



which implies that Eq. ( II 001) and Eq. (11051) are equal up to subleading terms. Indeed, 
if we integrate Eq. ( jl06l> we obtain 



01-w = a a l 



^ln 1 + (-x + 1) + i (-x 2 + l) + ...V 



(107) 
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This argument shows that at small x the 5 function in Eq. fllOOp cancels an anal- 
ogous contribution coming from the endpoint singularity of the plus prescription, 
thereby demonstrating that the first term in Eq. (11001) is not subleading in compar- 
ison to the second, and that small x leading terms in rapidity distributions may be 
unaccompanied by explicit lnx factors. In fact, there can be nontrivial cancellations 
among different terms due to endpoint singularity regularization. The situation is 
much clearer in Fourier-Mellin-space, where we can identify all small x enhanced 
terms as poles in N = or N = ±y. 

The emission from the other leg gives an identical contribution with only the 
replacement u — > 1/u, which in Fourier space amounts to b — > —b. The full two-leg 
NLO result in Fourier-Mellin space is thus 



du 



(N, b) = a a s 



+ 



N (N -~) N(N+%) 



(108) 



which in x-u space corresponds to the sum of Eq. fllOOp and the contribution from 
the other leg, obtained by letting in u — > 1/u in Eq. f II 1) . It is easy to check by 
inspection that this NLO result coincides with the small x limit of the full result 
in [25]. 

We now turn to the NNLO case, for which no closed-form analytic expression 
exists. We must compute the three diagrams of Fig. HI and recall that in the small x 
limit there are no interferences between them (see Sect. 12. 6j ). We start by considering 
the case of two emissions from the same leg, say from the upper leg (see Fig. @]-c). 
The result for the differential cross-section is El 
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l+e 



(109) 

The leading small x behavior in Eq. f 1 1 9 j) comes from the ordered region £ 2 £i, 
Zi 1. With this in mind we can perform the ^ integrations and the MS subtraction 
to obtain 



da\ 



2-gluons 
leg 



a a 2 s 



^ln 2 -5(ziz 2 - x) 
2 x 
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X {ZiZ 2 



X 



\yy{z x z 2 
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Z\Z 2 



X 



dz\dz 2} 
(110) 



3 From now on we will neglect all angular correlations, which can be shown to be subleading at 
small x. Consequently, as in the previous sections we will omit all angular terms like (47r) e /r(l — e), 
as they are always subtracted in the MS scheme. See Appendix \K\ for details. 
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which leads to the rapidity distribution 
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For the sake of comparison to the resummed result, it is convenient to consider 
the Fourier-Mellin transform of Eq. fillip . Starting from Eq. (I109p and performing 
the MS subtraction in Fourier-Mellin space we obtain 
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(112) 



Once again, despite naive appearance, the three terms in Eq. fillip are all of the 
same order. The small x expansion of Eq. f illip is 
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or 



j 2-gluons 
^1-leg 

dy 



x,y) 



8a a 2 yQlni-^ (l + 0(xe 2 ^)), (114) 

whose Fourier-Mellin transform is the leading-iV term given in the last step of 
Eq. dH2D. 

Having obtained the rapidity distribution for two emissions from the upper leg 
(Fig. 0]-c) it is immediate to obtain the result for two emissions from the lower one 
(Fig. HJ-a), through the substitution y — $> — y (i.e. u — > 1/u or b — > —b). To complete 
the small x NNLO distribution we only need the distribution with one emission from 
each leg. In this case the leading small x term in Fourier-Mellin space reads 



7 2-gluons 
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N 2 (N-ib/2) (N + ib/2)' 



(115) 



Adding everything up we get 

d<J NNLO 



dy 



(N,b)=a a 2 — 
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(116) 
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Combining the results of this section, we can write the small x rapidity distri- 
bution up to NNLO as: 



da , A T , . 



1 + 
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N \N -ib/2 N + ib/2 
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+ 0(a 3 s 



(117) 



Note that letting b = in this formula we reproduce the small x behavior of the 
inclusive cross-section. 



4.2 Rapidity distribution with finite m t0 p 

As already mentioned, the effective Lagrangian Eq. fl86l) is inadequate at small x 
and leads to spurious singularities. We now turn to the determination of the correct 
small x rapidity distribution. 

As before we start by considering the NLO distribution, i.e. we consider single 
gluon emission. For defmiteness, we consider first an emission from the upper leg. 
We obtain: 
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where 



(118) 



(119) 



and A\ is the form factor defined in [5] 0. Again, the hard part, i.e. the term in the 
first square bracket, reproduces the LO result when £ — > 0, z — > 1, as expected from 
Eq. ([15]): 



c^Lo(r) =a (127T 2 r) 2 |A 1 (0,0,r)| / = ( x (r). 



Also, because 
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12vr 2 



(120) 
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in the heavy top approximation Eq. (11181) reproduces our previous result Eq. ( 193|) . 

The form factor Ai acts as a cut-off for the £ integration at large £, thus removing 
the spurious small x singularities of the effective theory. To see how this works, we 
consider the rapidity distribution in Fourier-Mellin space. Starting from Eq. (11 IS)) 
we obtain: 
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(122) 



'Note that in Ref. [5] the definition of r is different. 
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where we have used the delta function to perform the N— Mellin transform. Without 
the form factor Ai (i.e. in the effective theory) the £ integral would be UV divergent 
at iV = 0. This leads to a spurious 1/N behavior of the integral. In the full theory 
this problem is no longer there: the integral is cut off by Ai, hence at small x we 
can safely let N = 0: 



d ^{N ) b) = ao(12vr 2 r) 2 
dy 



a. 



N - £6/2 



l+e 



(123) 



The £ integral now is well-behaved in the ultraviolet, but it still contains a 
collinear singularity. To extract it, we simply integrate by parts: 
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where we have used 
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since e < 0. We can thus write the MS subtracted rapidity distribution as 
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where we have defined 



|Ai(0,0,r)|' 



djM^r) 
d£ 



(127) 



Numerical results for ci(r) are tabulated in Tab. [TJ Note that the result Eq. (j!26p 
does not depend on N and b separately, but only on the combination N — ib/2. This 
implies that in (x, y) space the result has a rapidity delta function: 
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To obtain the result for emission from the lower leg it is sufficient to symmetrize 
this result. The full small x NLO result thus is 
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(129) 
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or in Fourier-Mellin space 



G^NLO 
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(130) 



Note that the small x behavior of the rapidity distribution with finite top mass 
is completely different from that in the effective theory: in the former case the 
rapidity is forced to sit at its endpoints, while in the latter we have a flat rapidity 
distributions (see Eq. (j!05p and recall dy = du/{2u)). Note also that the integration 
of Eq. (11291) leads to a constant, in contrast with the spurious lnx behavior of the 



effective theory. 

We can compare also the NNLO distributions. As in the previous subsection, 
we start from two emissions from the upper leg. In this case the small x rapidity 
distribution in Fourier-Mellin space is 
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(131) 

where we have regulated just one of the two collinear singularities coming from the 
extra gluon emission. To deal with Eq. ( I13ip we first note that we can safely let 
N = in the £ integral as before. To deal with the remaining collinear singularity, 
we just integrate by parts as in the NLO case. The MS result is then 
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(133) 



1^(0,0,^1^0 3 2 d£ 

Numerical results for c 2 are tabulated in Tab. [TJ 

The results with two emissions from the lower leg is trivially obtained by per- 
forming the replacement b — > — b in Eq. (I132p . In the nontrivial case of one emission 
from each leg the differential cross-section is 
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where A3 was defined in Ref. [5]. 

Note that if we let £ = in the first square bracket of Eq. (I134p and perform the 
angular integration we recover the one- leg result Eq. (11181) . From the differential 
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Table 1: Numerical values for the NLO coefficient C\ Eq. (I127P and for the NNLO 
coefficients c 2 Eq. (1 133[) and Eq. (11371) . 
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cross-section Eq. fjl34j) we can easily write down the rapidity distribution in Fourier- 
Mellin space 
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As in the previous case we can let iV = in the integral and perform the MS 
subtraction by integrating by parts. The final subtracted result now reads 
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Numerical values for c^i are tabulated in Tab. [TJ 

Adding up all contribution up to NNLO we obtain the following small x rapidity 
distribution: 
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5 Resummation of the Higgs rapidity distribution 

In this section we use Eq. ( IHTT) (or, equivalently, Eq. ( |84l) ) to perform the small x 
resummation of the Higgs rapidity distribution, both in the effective theory and for 
finite top mass. We will check that that up to NNLO the resummed results agree 
with the explicit computations presented in the previous section. 



5.1 Pointlike effective interaction 

The resummation of the inclusive cross-section in the limit m top — > oo was first 
performed in Ref. [I]. The resummed result for the gluon-gluon sub-process in 
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Mellin space is 
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where 7 = 7 (^f) is the LLx anomalous dimension. Because of the simple kinematic 
relation fl92|) between z and y, the only further piece of information needed for 
the resummation is the off-shell cross-section with the contributions from emissions 
coming from the two legs kept separate. This is given by |4J 
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where the two variables M, M correspond to radiation from either leg. Using 
Eq. (1841) . we immediately get 
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Note that, in accordance with Eq. (1841) . the argument N in the anomalous dimension 
is shifted by ±ib, but the intrinsic N dependence of the coefficient function is not. 

Recalling that -R^g = 0(af) it is immediate to check that the resummed result 
agrees with the fixed-order results up to NNLO obtained in Sect. 14.11 Eq. (I117p . 
Expanding Eq. (j!4ip up to NNLO and performing the inverse Fourier- Mellin trans- 
form, the result in terms of the partonic rapidity y is 
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5.2 Finite top mass 

We now consider the case of finite top mass. The resummation of the inclusive 
cross-section has been performed in Ref. [5]. The result was written in terms of a 
double Mellin transform 
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5 It follows that these results cannot be used to check a recent claim [35] that the standard result 
of R from Refs. [U [19] is incorrect. 
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where as above M and M correspond to radiation from either of the two legs. 
Expanding Eq. (j!43j) in powers of M and M one obtains the coefficient of the LLx 
singularity to any desired order in terms of double integrals over £ and £, which have 
to be evaluated numerically. Note that now in Eq. ( 11431) all the N dependence comes 
from j(a s /N), since in the small x limit the form factors Ai are iV— independent as 
one expects for a non-pointlike interaction. 

It is immediate to obtain the rapidity distribution: 
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where we have defined 
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Integrating Eq. (j!44p by parts we can make its perturbative expansion explicit: 
da n 



' 99 



dy 



(N,b) = a (r)R m (M b )R m (M b ) x 

{Ml + Ml) Ci(r) + Ml M " c ^) 

j,k>0 



i>0 
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(j-l)\(k-l)\ 

Integrating Eq. (I146p we obtain the right single logarithmic behavior of the total 
cross-section. 

For comparison with the explicit results of the previous section, we write the 
(x, y) resummed rapidity distribution up to NNLO: 



da gg (5(1 — x) 

—^(x,y) = <t (t) ' 

+a s c 1 (r) 



5[y--\a.x \ + 5 (y + ^lnx 
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which agrees with our previous findings. The pattern of Eq. (j!48p persists at higher 
orders: there are two endpoint contributions coming from emissions from just one 
leg, plus a bulk term coming from emissions from both legs. 

So far we have only considered the gluon-gluon partonic sub-process. However, 
beyond LO other channels open up. As discussed in Sect. 12.61 the high energy 
behavior of the quark initiated sub-processes can be straightforwardly computed 
from the gluon-gluon case thanks to color-charge relations. The LLx behavior for 
the inclusive Higgs production in the different partonic channels have been explicitly 
computed to NNLO in [9]. The generalization to the rapidity distributions is not 
difficult. Up to NNLO we have 



d [ y - -Inx) + S (y + -Inx 



ci(r 
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da, 



HI] in) 
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(x,y) = <t (t)< a l 



Cl>1 (r) + 0(a 3 s ) 



(149) 
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5.3 Matching to the effective theory 

As mentioned in the introduction, resummed results are useful not only for the sake 
of all-order resummed phenomenology, but also as a way of obtaining information 
on higher order perturbative terms. This is particularly interesting in the case of 
Higgs production, where results are most easily obtained in the m t0 p — > oo limit, 
which however fails at small x. It is then possible to improve the m top — > oo result 
by correcting its small x behavior with its exact finite m top form extracted from 
resummation. For the inclusive NNLO cross-section this was done in Ref. [6]. A 
more ambitious task is to construct an accurate approximation to the full finite 
m top result by matching the small x behavior from resummation to an expansion in 
powers of l/m top , which at the inclusive NNLO level was done in Ref. [§]. 

Here we perform a matched determination akin to that of Ref. [6] for the NLO 
rapidity distribution, by combining the small x behavior Eqs. (11291 11491) with the 
analytic expression in the m top — > oo limit from Ref. [25]. This quantity has also 
been obtained by numerical integration of a fully differential expression in Ref. [27J, 
with which we will compare our results. In Ref. [25] the rapidity distribution is 
parametrized by the variables z and y where z coincides with the variable that we 
have called x throughout this paper (compare in particular Eq. (|68|) ) while y should 
not be confused with the rapidity, and coincides with the variable w defined as 

w = 1 LZ£ (151) 
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in terms of x and u Eq. (|96|) . 

The matched result is constructed as follows: 



dcTij I 

dw leff 
da. 



ij_ I I dajj 

dw leff dw [matching 



X > X ma tch 
% — ^match 



(152) 



where ij=gg or qg, by ff we denote the result in the the heavy top limit as given 
in Eqs. (24)-(25) of Ref. [25], and -4^-1 , , . is a matching term which subtracts 

n v ' v ' 1 " dw I matching 

the spurious double log small x behavior from the heavy top result, and it replaces 
it with the correct small x behavior for finite mt op - 

The matching term in turn is determined by noting that at NLO at the inclusive 
level the spurious double logs correspond to terms which in Mellin space behave as 
either a simple or a double N pole (i.e. as functions of x which are either constant 
or grow as lnx as x — > 0). Thus, the matching terms are 



da, 



Ml 



dw 



da, 



qy 



dw 



matching 



matching 



X 



[w(x — 1) + 1] [w(x — 1) — x\ 
+ 3 Cl (r)[5(l-w)+6{w)}, 
x 2 -l 

[w(x — 1) + 1] [w(x — 1) — x] 
+ 2 Cl (r)[5{l-w)+6(w)}. 



- 3(2 - w(l - u>)) 



+ 



+ 

(153) 
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where the first line of Eqs. (I153|ll54p subtracts all contributions to the heavy top 
results (respectively Eq. (25) and Eq. (24) of Ref. [25]) which upon integration lead 
to contributions to the inclusive result which either grow or go to a constant at small 
x. The second line of both expressions adds back the correct single-log behavior as 
given respectively in Eq. (11291) . and Eq. (11491) . Note that the first term on the 
right-hand side of Eqs. (I153[ 11541) is the leading double-log contribution computed 
in Sect. 15. 1[ Eq. (I106j) . symmetrized, expressed in terms of w and multiplied by an 
appropriate Jacobian; when expressed in terms of the rapidity y this contribution is 
a flat rapidity distribution, as already noticed in Eq. ( 11301) . The second subtraction 
term on the right-hand side of Eqs. (I153[ I154j) is NLLx hence it is not determined by 
LLx resummation and we have extracted it from Ref. [25] directly. We will choose 
for a; mat ch the same matching value as used in Refs. [3 [6]. This ensures that our 
result reproduces the inclusive one once it has been integrated over the rapidity 
range. 

We have computed the full hadronic rapidity distribution up to NLO, using for 
the NLO term the matched results Eqs. (115211154]) . The hadron-level result is given 
by 



— = f dw f dx——C(x,w,Y), 
dY J J Xh dw 



(155) 



40 



where the partonic luminosity C is defined as 



Xh , / \/Xhe Y \ , ( JxJ^e 5 



C(x,w,Y) ee -£/ y-£= f , (156) 

and u is obtained by inverting Eq. ( I15ip . Results for two different values of the Higgs 
mass (m h = 130 GeV and m h = 280 GeV), at the LHC with VS = 7 TeV and VS = 
14 TeV are shown in Fig. [5] we have used NNPDF2.0 parton distributions [36] (PDF 
uncertainties are not shown). We plot the ratio of the NLO matched result to the 
large top mass result of Ref. [23], divided by the corresponding ratio of total cross- 
sections, so that only shape differences are shown in the plot. These plots confirm the 
conclusion of Ref. [27] that corrections to the NLO rapidity distribution due to finite- 
mass effects are below 5%. We also determine the precise shape of the correction, 
which could not be determined in Ref [27] because of insufficient numerical accuracy. 
In comparison to that reference, we also seem to find a somewhat more noticeable 
correction in the central rapidity region, though of course the effect on the total 
cross-section remains quite small. For heavier Higgs masses the correction to the 
shape becomes smaller. 

More refined matching procedures might be constructed, for example by taking 
the value of a; matc h to depend on y. However, the subtraction term (first line on 
the right-hand side of Eqs. (11531 11541) turns out to have a very small impact on the 
shape correction, at least for the values of the Higgs mass and energy considered 
here. This can be understood in part as a consequence of the fact that the leading 
double logarithmic term (first term on the right-hand side of Eqs. (I153[ I154j) is flat 
in rapidity, as mentioned above, so it does not affect the shape. 

A similar matched procedure cannot be carried out to NNLO because an analytic 
expression of the Higgs rapidity distribution to NNLO does not exist even in the 
heavy top limit, to the best of our knowledge. It is possible nevertheless to provide 
a rough estimate (more likely an upper bound) of the impact of finite top mass 
corrections on the NNLO cross-section by computing the i^-factor 

/da L0 da NL ° d<r™ LO \ I (da L ° Ar NLO da^ \ 

K ~\dY + dY + dY ) M / [dY + dY + dY J elI (157) 

where at NNLO level only the contribution from the small x tail is included, in turn 
approximated in each case by its LLx behavior, namely 

LLx - 1 dw dx^^C(x,w,Y) (158) 



dY J J Xh dw 

with dajj^j: /dw given by the NNLO term in Eq. (11481) for the matched case and 
by the NNLO term of Eq. (11421) for the effective (heavy top) case, in each case 
multiplied by the appropriate Jacobian. We found that at the LHC at 7 TeV the 
finite mass corrections are well below 1%, while at 14 TeV they are larger but still 
below 2% in the full rapidity range. This suggests that the use of the effective theory 
is fully justified for rapidity distributions also beyond NLO. 
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Figure 5: The ratio of the rapidity distribution for Higgs production computed up 
to NLO using the matched expressions Eqs. (I152H154]) to the result in the heavy top 
limit of Ref. [25]. The ratio is rescaled by the corresponding ratio of total cross- 
sections. Results are shown for the LHC at 7 TeV (left) and 14 TeV (right); the 
dashed yellow curves are for irth = 130 GeV and the solid green ones for rrih = 
280 GeV. 

6 Conclusions and Outlook 

In this paper we have accomplished for the first time small x resummation of a 
differential cross-section (as opposed to a fully inclusive one). This has been done 
by extending to differential rapidity distributions the so-called high energy or 
factorization, which expresses LLx cross-sections in terms of standard partonic cross- 
sections, but computed with incoming off-shell gluons. The result has been made 
possible by exploiting the duality which relates LLx and LLQ 2 evolution equations, 
to re-express high energy factorization in terms of standard collinear factorization, 
which in turns is straightforwardly extended to rapidity distributions. 

Our final resummation prescription is quite close to the well-established one 
derived long ago for inclusive cross-sections pQ . Indeed, it is based on evaluating the 
LO rapidity distribution for the relevant process but with incoming off-shell gluons. 
As in the inclusive case, the leading high energy behavior is found by means of a pole 
approximation in Mellin space, whereby the Mellin variable which is conjugate to 
the gluon virtuality is identified with the LLx anomalous dimension. However, one 
must now also take a Fourier transform with respect to rapidity, and the argument 
of LLx anomalous dimensions undergoes an impact-parameter dependent shift in 
the complex plane of the variable which is Mellin conjugate to the scaling variable 
x. 

As first application, we have performed a resummation of the rapidity distribu- 
tion for Higgs production in gluon-gluon fusion. We have computed the small x 
limit of the rapidity distribution to NNLO, both in the heavy top mass approxima- 
tion and for finite m top (only the m top — > oo NLO result being already available in 
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closed form). We have shown the way these expressions are related to the results 
obtained by expanding to finite order the resummation formula and inverting the 
Fourier-Mellin transform, and we have seen that this involves nontrivial cancellations 
between terms which at the resummed level have different origins, thereby providing 
a nontrivial consistency check on the resummation procedure. As a first example 
of application we have constructed an approximate analytic expression for the NLO 
rapidity distribution with finite top mass, which appears to be in agreement with 
the numerical results of Ref. [27], but provides a prediction which is not affected by 
problems of numerical accuracy. 

The potentially more interesting application of this formalism is to Drell-Yan 
rapidity distributions, which will be explored at very small x at the LHC, specifically 
by the LHCb collaboration [18], and which may play an important role both in 
determining parton distributions and testing resummed QCD predictions at small 
x. More in general, our result provide a first generalization of small x factorization 
and resummation beyond the inclusive LLx level. It will be interesting to see whether 
they can be generalized not only to less inclusive quantities, but also beyond the 
leading logarithmic approximation. 
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Appendix 



A The MS subtraction procedure to all orders 

Here we briefly sketch how to subtract collinear singularities arising from kernel 
iteration within the MS scheme. Here we outline only the features which are relevant 
to understand how Eqs. ( l3"Tj) . (j4"2"j) are obtained. For a general treatment, including 
a discussion of ultraviolet and infrared singularities, as well as for a more detailed 
explanation, we refer the reader to the original papers [301 121] ■ 

We start by writing a generic cross-section as the convolution product of a "bare" 
(i.e. unsubtracted) coefficient function M and a "bare" PDF f: 

a = M(g)f, (159) 

where <g> denotes the standard convolution product in x space. We assume to work 
in d = 4 — 2e dimensions in order to regularize M and f. Following Sect. |2j we now 
write the coefficient function M as the product of a hard part H and the iteration 
of a kernel K: 

M = H® x>kT {l+K+K® x>kT K+...) = M(l+K+K 2 +K 3 + ...) = H——, (160) 

1 — K 

where now the convolution ® Xy k T stands both for normal convolution in x space 
and for kx integration in the transverse space. This convolution product will be 
understood in the following formal manipulations. 

As a result of the kx integrations, the coefficient function Eq. f)160p contains 
multiple collinear poles of the form 1 / e J and thus it is meaningless in 4 dimensions. 
To perform the subtraction, we first introduce the "pole-part" projector P which 
acting on a generic function F selects just its l/e k poles, i.e. 

fc>0 

We now use this projector to factorize the kernel 1/(1 — K) in Eq. (I160p into a finite 
part and a pure-pole part. To this purpose, we first consider the kernel (1 — K) and 
perform the following formal manipulations: 

1-K = 1-PK-(1- P)K = [1 - PK(1 - (1 - P)K)- 1 ] [1 - (1 - P)K) , (162) 

which leads to (note the reverse order): 



1 



K 



1 



1 - (1 - P)K 



l-PK{\ - (l-P)K)- 1 



(163) 



This is the desired factorized expression: the term in the first square bracket is 
by construction free of collinear singularities, which have been moved to the second 
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term. Using this factorized form, we can now rewrite the cross-section f 1 1 5 9 j) as 



a = MT 
with 



H 



1 - (l-P)K 



1 - PK{1 - (1 - P)K)~ 1 



r 



Mr, 



(164) 
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r 



H 



1 - (l-P)K 
1 



l-PK{\- (1 -P)K)- 1 



r. 



(165) 



Now both M and T are finite, hence we can safely let e = and obtain a well-defined 
d = 4 result. Note that this subtraction scheme is uniquely defined once the action 
of P is specified. The subtraction of the pure pole part (as in Eq. (11611) ) defines the 
MS scheme. In the MS scheme one chooses to subtract also terms proportional to 
—je + ln47r coming from angular integrations. At leading logarithmic accuracy the 
action of P is defined as: 



Pm F = ekp ) x ]k ex P [* e (-7b + ln 

k>0 



(166) 



We now illustrate this subtraction procedure on the simple case of the t— channel 
iteration of Altarelli-Parisi gluons. In this case we have 
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and so on. Note that K 2 contains both single and double poles: 



K 2 

"ho 



! In (£) + ln(4vr) - 7s 



2e 2 
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x ln 
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(168) 



Careless subtraction of all poles and ln(47r), 7e terms would lead to the incorrect 
result 

s 2 u b = « 2 7o(A01n 2 ^. (169) 
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The correct MS result is instead obtained following the procedure explained above: 
instead of subtracting poles from K 2 we should consider (1 — P)[K(1 — P)K\. In 
this case we get 

( 1 (4tt) 
(l-P)K = (a sl0 (N))l 

K{l-P)K = (a sl0 (N)) 2 I t^T^T— ( ^ ) + < -™ 



P) [K(l - P)K] 




K7o(iV)r ^ ( 7* ) 1 : 'I'D 



with S € = exp [e (ln(47r) — 7s)]- Note that this time 

1 - P) [K{1 - P)K] = i- (a sl oY In 2 ^ + 0(e), ( I 72 ) 



^(« s7 o) 2 ln 2 ^ 



which is the correct MS result. 

This subtraction procedure can be easily iterated. In particular, with n kernels 
the result is 

= (a sl0 r-\n n % + O(e). (173) 

Note that before performing the last (1 — P) subtraction, the n-kernel result contains 
only a single l/e n pole. 

A particularly important feature of this subtraction method is its iterative na- 
ture. This is especially well suited for the purpose of this paper, because we can 
perform all the universal kernel subtractions once for all, and consider separately 
the single collinear divergence coming from attaching the universal ladder to the 
process-dependent hard coefficient function, see Eqs. (I3T|) . f l4"2|) . Also, the MS sub- 



traction factor S t exactly cancels all terms coming from (47r) e /r(l — e): the full MS 
result can be obtained by just neglecting all (47r) <: /T(l — e) terms and using S t = 1. 
This simplification only relies on the fact that each emission is accompanied exactly 
by a (47r) e /r(l — e) term, which happens if there are no angular correlations. Since 
in all our derivations angular correlations are subleading, in the main paper we have 
used this simplified version of the subtraction procedure. 
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B Resummation of rapidity distributions in MS 



In this appendix we provide an explicit derivation of the resummation of the ra- 
pidity distribution in the MS scheme. Because rapidity is determined by the lon- 
gitudinal momentum components, while collinear singularities appear in transverse 
momentum integrations, the main features of the MS derivation are the same of the 
derivation in Sec. |3j 

For simplicity we consider emissions from just one leg. Extension to the full case 
is straightforward. We start by writing the n— kernel fully differential cross-section 
in d = 4 — 2e dimensions, see Eq. (|4ip |j: 

d(T Q^' ^' 

p ' as] e ) ir^fe x c z *> as] e ) x 

(174) 

where as in Eq. (!73|) the splitting function P is the inverse Mellin transform of the 
anomalous dimension 7 in Eq. ( 14"T|) . Note that in Eq. f |174j) we have written C as the 
most generic (2PI) function of all the relevant variables. Quasi-collinear kinematics 
determines the dependence of C on its arguments. Because rapidity distributions do 
not factorize in N— Mellin space, for the time being we write everything in x— space. 

The d— dimensional rapidity distribution is immediately obtained from the fully 
differential cross-section Eq. (II 74ft : 

da n ( fi 2 
~dy V ' Q 2 ~ ,y,as,e 



x 



f°° d ^n , c \ f d ^~i f d ^ 

/ 7ir- e Cy(x,Zi,^,y,a s ;e) __ x ... x / -j 

JO Sn JO ? n -l JO SI 



l+e> 



(175) 



where the rapidity-dependent coefficient function is defined as 

C y (x,z i ,Z i ,y,a.ie) = --± / — / dU T \M (V(n) + g*(p L ) -> ^)| 2 x 

x5 4 (Pj -Pf)x8(y-\\n ^±^1) , (176) 



6 Here again we omit all (47r) e /r(l — e) terms in view of the MS subtraction, as discussed in 
Appendix |A] 
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in terms of the spin and color averaged squared amplitude \M.\ 2 for the off-shell 
process V(n) + g*(pi) — > J 7 (see Fig. |2J), and the polarization sum for the off-shell 
gluon is performed through the projector V Eq. (Tl2|) . We now use the fact that we 
are in quasi-collinear kinematics (i.e. £i <C £»+i) an d in the small x regime Eq. flSJ). 
This implies first, that C y is sensitive only to the largest transverse momenta £ n 
and second, that up to terms which vanish when k\ <C s the effect of the ladder 
insertion to C y is just a longitudinal boost. We have thus 

C v (x,Zi,£i,y,at s ;e) = 
x 



>Zn,y- o ln ^i 

Z\...Z n I 



-lnz n ,a s ;e )+0 ^ . (177) 



Since upon integration 0(k^/s) terms give rise to contributions suppressed by pow- 
ers of x, we can safely discard them. 

Thanks to quasi-collinear kinematics we can thus write Eq. f )175p as 
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(178) 



From now on we will systematically omit all 0(x) terms. 

Convolution products in Eq. (I178P turn into ordinary ones by Fourier-Mellin 
transformation: 
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(179) 



This expression is identical to Eq. f l4"Tj) . but with the anomalous dimension 7 is 
evaluated at N — ib/2. This does not interfere with collinear singularities, hence the 
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MS subtraction is performed as in the inclusive case, with the result (see Eq. (|44|) ) 
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where 7^ was defined in Eq. (14"5|) . As in the inclusive case, Eq. (I180p has a finite 
e — 7- limit: 
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Using Eq. (j!8ip the resummed result immediately follows (with fi 2 = Q 2 , thus 
omitting the explicit dependence on fi 2 /Q 2 ): 



da 
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which is the same of Eq.t HTl) but this time with a shifted 1Z factor. Note that at 
LLx we also have 
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We finally show how to deal with a scheme change in rapidity distribution, 
thereby in particular explaining the shift in the argument of 1Z in Eq. (I182p . To this 
purpose, assume a scheme change from MS to some "small 2" scheme is performed 
by dividing the parton distribution by a LLx function M\ 



Msmall x = MygrJ\f ) 
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Using the MS kernel expansion Eq. (1165p we can rewrite M sma u x as 
1 



small x 
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\-{l-P)K 



N = H [1 + (1 - P)K + {1-P) [K{1 - P)K] + ...] M. 

(185) 
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where (see Appendix [A} all products must be understood as convolution in x and 
kx spaces. It follows that the effect of the scheme change is the same as that of a 
kernel insertion on the first (outer) rung of the ladder: 
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In Fourier-Mellin space the convolutions in Eq. (11861) turn into ordinary products 
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so the argument of Af ends up being also shifted when used for computing rapidity 
distributions. Clearly this holds also if we consider the full ladder and not only a 
single kernel insertion. 

Collecting the results of this section, we can write the full MS rapidity distribu- 
tion resummation formula as 
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(188) 



The extension of Eq. (11881) to emissions from both legs along the lines of Sec. [3] is 
straightforward 
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(189) 



50 



References 

[1] S. Catani, M. Ciafaloni and F. Hautmann, Nucl. Phys. B 366 (1991) 135. 

[2] S. Catani and F. Hautmann, Nucl. Phys. B 427 (1994) 475 [arXiv:hep- 
ph/9405388]. 

[3] R. D. Ball and R. K. Ellis, JHEP 0105 (2001) 053 [arXiv:hep-ph/0101199]. 

[4] F. Hautmann, Phys. Lett. B 535, 159 (2002) [arXiv:hep-ph/0203140]. 

[5] S. Marzani, R. D. Ball, V. Del Duca, S. Forte and A. Vicini, Nucl. Phys. B 
800, 127 (2008) [arXiv:0801.2544 [hep-ph]]. 

[6] S. Marzani, R. D. Ball, V. Del Duca, S. Forte and A. Vicini, Nucl. Phys. Proc. 
Suppl. 186 (2009) 98 [arXiv:0809.4934 [hep-ph]]. 

[7] S. Marzani and R. D. Ball, Nucl. Phys. B 814, 246 (2009) [arXiv:0812.3602 
[hep-ph]]. 

[8] G. Diana, Nucl. Phys. B 824, 154 (2010) [arXiv:0906.4159 [hep-ph]]. 

[9] R. V. Harlander, H. Mantler, S. Marzani and K. J. Ozeren, arXiv:0912.2104 
[hep-ph] . 

[10] R. V. Harlander and K. J. Ozeren, JHEP 0911 (2009) 088 [arXiv:0909.3420 
[hep-ph]]. 

[11] A. Pak, M. Rogal and M. Steinhauser, JHEP 1002 (2010) 025 [arXiv:091 1.4662 
[hep-ph]]. 

[12] L. N. Lipatov, Sov. J. Nucl. Phys. 23 (1976) 338 [Yad. Fiz. 23 (1976) 642]. 

E. A. Kuraev, L. N. Lipatov and V. S. Fadin, Sov. Phys. JETP 45, 199 (1977) 
[Zh. Eksp. Teor. Fiz. 72, 377 (1977)]. 

I. I. Balitsky and L. N. Lipatov, Sov. J. Nucl. Phys. 28, 822 (1978) [Yad. Fiz. 
28, 1597 (1978)]. 

[13] V. S. Fadin and L. N. Lipatov, Phys. Lett. B 429 (1998) 127 [arXiv:hep- 
ph/9802290]. 

G. Camici and M. Ciafaloni, Phys. Lett. B 412 (1997) 396 [Erratum-ibid. B 
417 (1998) 390] [arXiv:hep-ph/9707390]. 

[14] See M. Ciafaloni, D. Colferai, G. P. Salam and A. M. Stasto, Phys. Rev. D 68 
(2003) 114003 [arXiv:hep-ph/0307188] and ref. therein. 

[15] See G. Altarelli, R. D. Ball and S. Forte, Nucl. Phys. B 742 (2006) 1 [arXiv:hep- 
ph/0512237] and ref. therein. 



51 



[16] R. D. Ball, Nucl. Phys. B 796 (2008) 137 [arXiv:0708.1277 [hep-ph]]. 

[17] G. Altarelli, R. D. Ball and S. Forte, Nucl. Phys. B 799 (2008) 199 
[arXiv:0802.0032 [hep-ph]]. 

[18] R. McNulty, In *Perugia 2008, Multiple Partonic Interactions at the LHC 
(MPI08)* 155-163 

[19] S. Catani, M. Ciafaloni and F. Hautmann, Phys. Lett. B 307 (1993) 147. 

[20] G. Altarelli, R. D. Ball and S. Forte, Nucl. Phys. B 575, 313 (2000). 

[21] G. Curci, W. Furmanski and R. Petronzio, Nucl. Phys. B 175 (1980) 27. 

[22] R. D. Ball and S. Forte, Phys. Lett. B 359 (1995) 362 [arXiv:hep-ph/9507321]. 

[23] M. Ciafaloni, Phys. Lett. B 356 (1995) 74 [arXiv:hep-ph/9507307]. 

[24] G. Altarelli, R. D. Ball and S. Forte, Nucl. Phys. B 599 (2001) 383 [arXiv:hep- 
ph/0011270]. 

[25] C. Anastasiou, L. J. Dixon and K. Melnikov, Nucl. Phys. Proc. Suppl. 116 
(2003) 193 [arXiv:hep-ph/0211141]. 

[26] C. Anastasiou, K. Melnikov and F. Petriello, Nucl. Phys. B 724 (2005) 197 
[arXiv:hep-ph/0501130]. 

[27] C. Anastasiou, S. Bucherer and Z. Kunszt, JHEP 0910 (2009) 068 
[arXiv:0907.2362 [hep-ph]]. 

[28] T. Jaroszewicz, Phys. Lett. B 116, 291 (1982); 

[29] M. Ciafaloni, Nucl. Phys. B 296 (1988) 49. 

[30] R. K. Ellis, H. Georgi, M. Machacek, H. D. Politzer and G. G. Ross, Nucl. 
Phys. B 152 (1979) 285. 

[31] M. Ciafaloni and D. Colferai, JHEP 0509 (2005) 069 [arXiv:hep-ph/0507106]. 

[32] S. Marzani, R. D. Ball, P. Falgari and S. Forte, Nucl. Phys. B 783 (2007) 143 
[arXiv:0704.2404 [hep-ph]]. 

[33] M. Kramer, E. Laenen and M. Spira, Nucl. Phys. B 511 (1998) 523 [arXiv:hep- 
ph/9611272]. 

[34] C. Anastasiou, S. Bucherer and Z. Kunszt, JHEP 0910 (2009) 068 
[arXiv:0907.2362 [hep-ph]]. 

[35] R. Kirschner and M. Segond, Eur. Phys. J. C 68 (2010) 425 [arXiv:0910.5443 
[hep-ph]]. 

52 



[36] R. D. Ball, L. Del Debbio, S. Forte, A. Guffanti, J. I. Latorre, J. Rojo and 
M. Ubiali, Nucl. Phys. B 838 (2010) 136 [arXiv: 1002.4407 [hep-ph]]. 



53 



